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Abstract 



We construct a radiative inverse seesaw model with local B — L symmetry, and investigate 
Ql| the flavor structure of the lepton sector and the fermionic Dark Matter. Neutrino masses 

^-i are radiatively generated through a kind of inverse seesaw framework. The PMNS matrix is 

— ■ derived from each mixing matrix of the neutrino and charged lepton sector with large Dirac 

CP phase. We show that the annihilation processes via the interactions with Higgses which are 
\ independent on the lepton flavor violation, have to be dominant in order to satisfy the observed 

\ relic abundance by WMAP. The new interactions with Higgses allow us to be consistent with 



the direct detection result reported by XENON100, and it is possible to verify the model by the 



|^ ■ exposure of XENON100 (2012). 
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1 Introduction 



Inverse seesaw mechanism which generates neutrino masses due to small lepton number violation is 
one of the intriguing way to describe tiny neutrino masses [U[2]. Thus interesting phenomenologi- 
cal implications have been accommodated [3H5]. However Dark Matter (DM) candidate has to be 
introduced independently if we discuss DM phenomenology in this kind of models. On the other 
hand, radiative seesaw mechanism is possible to relate neutrino mass generation with the existence 
of DM [6HH] . In particular, the radiative seesaw model which is proposed by Ma [6] is the simplest 
model with DM candidates. Subsequently there are a lot of recent works of the model [SHU] and 
the extended models [T2HT9] . The other radiative neutrino mass models are studied in Refs. [2DT423] . 
The Z2 parity imposed to the model forbids to have the Dirac neutrino masses, produces neutrino 
masses at one loop level, and stabilizes DM candidates. Therefore the feature of the radiative seesaw 
model motivates to connect the existence of DM with the neutrino mass generation due to inverse 
seesaw mechanism. 

In this paper, we construct a radiative inverse seesaw model with U(1)b-l as a concrete example 
and analyze the neutrino masses, mixing and the feature of DM. We add three pairs of B — L charged 
fermions and a scalar to Ma model [5] . The scalar particle breaks the B — L symmetry spontaneously. 
In the model, the small Majorana mass terms which violate U(1)b-l weakly and explain the tiny 
neutrino masses are generated from a higher operator when the U(1)b-l symmetry is spontaneously 
broken. We assume the structure of neutrino mass matrix that induces the best fit value of 8 12 of 
the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino mixing matrix [21]: sin 2 # 12 =0.311, which 
is within la confidence level in the global analysis [25]. Moreover, we introduce the charged-lepton 
mixing (A) with nonzero Dirac CP phase to induce non-zero #13 recently reported by T2K [26], Double 
Chooz [27], Daya-Bay [28], and RENO [29]. This method was firstly proposed by S. King [30~|[3"T|. 
Due to the mixing A, neutrino Dirac Yukawa couplings are strongly constrained by the lepton flavor 
violation (LFV), especially, \i — > e-f process. Hence in the original radiative seesaw model [6], it 
is hard to derive the observed relic density of DM [32] associated to the annihilation channel if we 
assume that the lightest right-handed neutrino is DM. However since the radiative inverse seesaw 
model has the other channels with different coupling due to the additional Higgs boson coming from 
the B — L symmetry, the observed relic abundance can be naturally obtained via the channel. Since 
the additional Higgs boson mixes with the SM like Higgs, the direct detection with Higgs portal also 
comes into the target of a discussion whether the result can be consistent with the experimental limit 
reported by XENON 100 [33] • It implies that the model can be investigated in the direct search of 
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DM near future. Therefore it is favored to be Higgs portal DM in the radiative inverse seesaw model 
from the consistency with neutrino masses, mixing, LFV and the DM property. 

This paper is organized as follows. In Section 2, we construct the radiative inverse seesaw model 
and its Higgs sector. In Section 3, we discuss the constraints from LFV, especially, [i — > process 
and the DM relic abundance. In Section 4, we analyze the direct detection of our DM including all 
the other constraint. We summarize and conclude the paper in Section 5. 

2 The Radiative Inverse Seesaw Model 
2.1 Neutrino Mass and Mixing 
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Table 1: The field contents and the charges. Three generations of right handed neutrinos N? and 
additional fermions Si, S^ are introduced. A pair of fermions S and S' is required to cancel the 
anomaly. The SM, inert and B — L Higgs bosons are denoted by $, 77, and \i respectively. 

The radiative inverse seesaw mechanism can be realized by introducing the B — L symmetry 
which is spontaneously broken at TeV scale ^j. The field contents of our model is shown in Table 1 
and the Lagrangian in the lepton sector is 

C = y e L$e c + y„L v N c + y s N c X S + ^ 2 S 2 + V + h.c, (2.1) 

where A is a cut-off scale and the generation indices are abbreviated. Note that the mass term SS' 
can be fo rbl dde n by the Z 2 synnnetj We have another dnnension 5 operator s nch as L*S V , 
which affects on the neutrino mass, and it would be difficult to forbid it. Thus we need to assume 
the coupling of this operator is small enough. 

After the symmetry breaking, that is x = (x° + u ')/V^Ll = (0° + V )/V2 with $ = (</> + ,</>) T 

1 Another way to induce the mass term of S is proposed by E. Ma [M], in which the term is done at loop level. 

2 The value of v' can be constrained by the LEP experiment that tells us mz>/g' = \Y^_ L \v' > 6 TeV [55] . where 
mz' and g' are the B — L gauge boson and the B — L guage coupling, respectively. Since Yg_ L =-l/2 is taken in our 
case, v' > 12 TeV is obtained. 
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and v = 246 GeV, the neutrino sector in the flavor basis can be written as 



C u m = y u L v N c + M T N C S + + h.c, 



Z 1 c2 



(2.2) 



where M = ysv' /y/2 and /i = Xsv' 2 / A. The scale of /x ~ 1 keV which corresponds to A ~ 10 14 GeV 
is required as we will see in the following section. The inert doublet r] = (r] + , (t)r + ir]j)/\/2) T does 
not have any vacuum expectation values (VEV). As a result, the 6x6 neutrino mass matrix in the 
basis (N c , S) takes the form: 

M T 
M fi 

The mass matrix M and \i cannot be diagonalized simultaneously in general. We assume that the 
mass matrix M and /i are diagonal as M = diag (Mi, M 2 , M 3 ) and /i = diag (/ii, /i 2 , Hz) for simplicity. 
Then the 6x6 mass matrix }A V can be diagonalized for every family as diag(m i+ , mj_) by the unitary 
matrix Ui for i-th family. The mass eigenvalues are expressed as 



(2.3) 



m i± 



Mf + — ± — . 



4 2 

When /ij <C Mi, the neutrino masses are degenerated. The unitary matrix Ui can be expressed as 

Mi iMj 



(2.4) 



Ui 



/ M 2 +m 2 / M 2 +m 2 



(2.5) 



where (mj+,mi_) = Uf M. l u Ui and .MJ, implies the 2x2 mass matrix for i-th family. The flavor 
eigenstates are rewritten by the mass eigenstates Vi± as follows: 
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(2.7) 



The light neutrino mass matrix seen in Fig. [T] is given as Ref. [6] by 1-loop radiative correction: 
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Figure 1: Neutrino mass generation via radiative inverse seesaw. 



where rriR and rrij imply masses of t]r and m . When fa <C Mi, we can obtain the approximate light 
neutrino mass matrix by expanding mj± up to the leading order, 
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where we will define the function and the parameter Aj as 
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2(4vr) 2 

We can see from Eq. (12.91) that the mixing matrix of the light neutrino mass matrix is determined 
by the structure of the neutrino Yukawa matrix y v since the majorana mass matrix \x is assumed to 
be diagonal in this case. In order to identify the structure of y u , here we set a specific texture of the 
neutrino mass matrix, which induces the best fit values of #12, as 

(a b -b \ 

m u = B (3A + VSB)/Q -(3A + V3B)/6 . (2.11) 
^-B -(3A + V3B)/6 (3A + V3B)/6 j 

Then the mass matrix can be diagonalized by the following mixing matrix and the eigenvalues are 
written as 

< y/( l + l/V7) /2 -y'^V^ N 

'l - 1/V7/2 -Jl + l/y/7/2 1/V2 , (2-12) 
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\ y/l + l/V7/2 l/y/2 J 



mi 



m 2 = A + — = (l + V7^j , m 3 = 0. 



(2.13) 



4 



Thus the squared mass differences are Am s 2 ol 



rrin 



m 2 and Am 



2 

atm 



m\ — m 



31) 



and the neutrino 



mass hierarchy is predicted to be inverted. In order to get Am^ = 7.62 x 10 5 eV 2 , Am 2 tm 
2.40 x 10 -3 eV 2 , which are the best fit values [25], we find the following solutions: 



(A, B) = (±4.92 x 10 -2 , ±2.53 x 10" 



;±1.83 x 10" 2 , ^3.23 x 10~ 2 ) eV. (2.14) 



Notice that the other solutions do not exist any more. We cannot obtain non-zero sin #13 from the light 
neutrino mass matrix Eq. ( 12. lip . Non-zero sin # 13 is derived from the charged lepton mixing as we will 
discuss below. We find J2i m « — 7.7 x 10~ 4 eV, and the effective mass (m ee ) = | £\ (O v l] 



mi 



0.026 



eV and tan #12 = (l — \/7) /\/6 which gives sin 2 9 12 = 0.311 which is the best fit la value [25]. The 
recent experimental value for (m ee ) are referred in Ref. [36J. We can choose the following texture 
which leads the above neutrino mass matrix^ : 
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(2.15) 



where the parameters a, b, c are expressed by A and B as follows: 
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(2.16) 
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To induce non-zero 9 13 , we consider the charged lepton mixing [301131] . If we set the Dirac mass 
matrix of charged leptons m e and the mixing matrix U e i as 

/ iA 7/^ n \ / 1 \e i5 \ 

eL ~ -\e~ iS 1 

V l ) 

where we define (|m e | 2 , Ittt,^ | 2 , |m T | 2 ) = U\ L m e m\U e L and 5 is the Dirac CP phase. From the mixing 
matrix O v l and U e L, we can obtain each of the element of the PMNS matrix, which is defined as 
Upmns = U\ L O v P, is found by 

A 



vi vl 

V yi j 



sin ^13 ~ sin 9 12 

v2 



(1 - l/V7)/2 + A=\f( l + 1/^/2 cos5, sin# 23 ~ -^=, (2.18) 



3 Our paramctrization is taken so that DM v\_ couples to the charged leptons. Even if one selects another 
parametrization, e.g, by replacing the first column and the second one of y v , the severe constraints from LFV do 
not change as we will discuss in the next section. 
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Figure 2: sin #13 versus CP Dirac phase 5. Here we restrict sin 9 12 within the range of the best fit 
with la that is 0.303 < sin 2 # 12 < 0.335 [25] . 

where P contains two majorana phases. The allowed value of sin # 13 is shown in Fig. [2] as the function 
of 5, within the range of the best fit with la that is 0.303 < sin 2 6*12 < 0.335 [25] . The light red 
region is in la range of sin 9i2. It suggests that large CP Dirac phase is required to satisfy the current 
global experimental limit of sin #13 for inverted hierarchy such as 0.023 < sin 2 #13 < 0.030 [25] which 
is shown as the two green sandwiched regions. 

2.2 Higgs Sector 

The Higgses <f)° and x° mrx after the symmetry breaking and these are not mass eigenstates. Inter- 
actions should be written by mass eigenstates in order to analyze the DM relic density and direct 
detection in the next section. The Higgs potential of this model is given by 

V = m 2 $ f $ + mlrfr] + m 2 x f X + Ai($ f $) 2 + A 2 (r/ t ?7) 2 + A 3 ($ t $)(?7 t ?7) + A 4 ($ t ?7)(?7 t $) 

+A 5 [($^) 2 + h.c] + A 6 ( X f x) 2 + A 7 (x t x)($ t $) + Xsi^xWv), (2-19) 

where A5 has been chosen real without any loss of generality. The couplings Ax, A2 and Xq have to 
be positive to stabilize the potential. Inserting the tadpole conditions; m 2 = —Xiv 2 — XjV 12 /2 and 
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-Xqv' 2 — Xyv 2 /2, the resulting mass matrices are given by 
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(2.20) 
(2.21) 

(2.22) 

(2.23) 



where /i implies SM-like Higgs and H is an additional Higgs mass eigenstate. The tadpole condition 
for 77, which is given by ^ 



VEV 
2 



0, tells us that 



mi, > 0, A 2 > 0, A 3 + A 4 + 2A 5 > 0, A 8 > 0, 



(2.24) 



in order to satisfy the condition v v = at tree level. The masses of 0° and x° are rewritten in terms 
of the mass eigenstates of h and H as 



<p = hcosa + H since, 
X° = —hsina + Hcosa. 



(2.25) 



3 The Constraints from Lepton Flavor Violation and DM 
Relic Density 

3.1 Lepton Flavor Violation 



We investigate Lepton Flavor Violation (LFV) under the flavor structure Eq. (I2.15p . We put a 
reasonable approximation m>i± ~ M, hereafter, since the scale of [i that is keV scale is negligible 
compared to the scale of Mj that is expected O(l0 2 ~ 3 ) GeV. The experimental upper bounds of the 
branching ratios are Br (// — > ej) < 2.4 x 10~ 12 [37], Br (r — > ^7) < 4.4 x 10~ 8 and Br (r — »■ ej) < 
3.3 x 10~ 8 [38]. The branching ratios of the processes £ a — > £37 (£ a ,£p = e, r) are calculated as 
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7 



where a cm = 1/137, Br (/i — > eu e u^) = 1.0, Br (r — > eu e u T ) = 0.178, Br (r — > iiv^v^) = 0.174, M n is 
rj + mass, Gp is the Fermi constant and the loop function F 2 (x) is given by 

„ , . 1 - 6x + 3x 2 + 2x 3 - 6x 2 In x /n „, 

F2(X) = 6(1^ • (3 - 2) 

The \x — > e7 process gives the most stringent constraint. If the mixing matrix of the charged leptons 
Eq. (I2.17P is diagonal, i.e. A = 0, the \i — > process does not constrain the model since the 
branching ratio Eq. (13.11) can be zero. Although r — > and r — > &y processes remain as LFV 
constraint, these are much weaker than /i — > ej. Instead of that, non-zero 6* 13 is not derived. Thus 
we can insist that fi — > e'-f constraint is closely correlated with non-zero #13. In order to be consistent 
with LFV and obtaining non-zero #13, the neutrino Yukawa couplings must be small enough to escape 
the LFV constraint. 



3.2 DM Relic Density 

If we assume that DM is fermionic and mass hierarchy Mi < M 2 < M 3 for the right-handed neutrinos 
Ni, the mass eigenstates v\± can be highly degenerated DMs due to the weak lepton number violation 
term ^ which induces small neutrino masses. Thus we have to take into account coannihilation of 
vi- and v\ + . The typical interacting terms are found as 



C 



hit 



(U eL y u ) al (£ a 7] + - is a r]°) 
. (ZAs)nSina 



i 1 



2 

Df cos a 



-h (fi 2 _ + 



V2 Vz 

(ys) u cos a 



+ h.c. 

H (v\_ + 



<l+) 



hff 



Df sin a 



Hff, 



(3.3) 



y/2 V2 

where the masses of rf^ assumed to be always heavier than v± to avoid the too short lifetime of 
DMs through our analysis. Three types of the coannihilation processes exist and these are shown in 
Fig EJ The effective cross section to ££, ff and hh are given as 
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Figure 3: t, u and s-channel of coannihilation processes of DM i/i±. 



where rjR, rjj and n + masses are regarded as same parameter M v for simplicity and only top pair is 
taken into account in fermion pair // because of the largeness of the Yukawa coupling. The SM-like 
Higgs mass and decay width are fixed to m>h = 125 GeV and Th = 10~ 2 GeV, and the heavy Higgs 
mass is assumed to be tuh < 200 GeV. The decay width of the heavy Higgs Th is expressed as 

y 2 t sin 2 a ( 4m 2 \ 3/2 
T// = — m H 1 r . 3.7 

The contribution of the process H — > vi±i>i± is also added to the decay width when the relation 
rrin > 2Mi is satisfied. Note that the mixing matrix of charged leptons U e L is introduced in the 
effective cross section a eS v Te \ since the initial Lagrangian (12. ip is not assumed as diagonal base of 
charged leptons. The s-wave vanishes and p-wave only remains in the above annihilation cross section 
due to the helicity suppression. 

Since the neutrino Yukawa couplings y v are severely restricted by LFV, the annihilation cross 
section of ££-channel is too small to obtain the proper DM relic abundance. Thus large contributions 
from tt and hh channels are required to have sizable effective annihilation cross section. These 
processes are possible because of the mixing of Higgses 0° and x°> namely the DMs v\± can be Higgs 
portal DMs. This is a different aspect from the loop induced neutrino mass model [6] and the various 
analysis of the right-handed neutrino DM in the original model [9| [39H42j . 

The independent parameters which appear in the analysis are A 1; A 3 , M 1; M 3 , M v , m H , (ys)n 
and sin a. The parameters A and B of the neutrino mass matrix Eq. (12. lip are determined by 
neutrino mass eigenvalues as Eq. (I2.14p . and we take the fourth solution of Eq. (12.141) as example. 
The parameter A of the charged lepton mixing matrix is fixed by the experimental value of sin 6*13 
as Eq. (I2.18p . and A 2 and M 2 are not relative with the analysis since the second column of the 
neutrino Yukawa matrix Eq. (I2.15P is zero. The allowed parameter spaces from LFV and the DM 
relic abundance for {ys)u — 0.5 and 1.0 in Mi-rriH plane are shown in Fig. |4] where the parameters 
are fixed to A, = 1 eV, M 2 = 1.5 TeV and M3 = 2.0 TeV. The parameter choice A« ~ 1 eV means 
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Figure 4: The allowed parameter spaces for (ys)u — 0.5 and 1.0 in M\-mn plane. The other 
parameters are fixed to Aj = 1.0 eV, M 2 = 1.5 TeV, M 3 = 2.0 TeV. The unfixed parameters are Mi, 
M v , mu and sina. 



that Hi ~ 1 keV if / (Mf/m 2 R ) ~ / (Mf/rrij) ~ 0.1 is assumed as can be seen from Eq. O2.10p . The 
result does not practically depend on Aj, M 2 and M3 since the dependence of Aj appears in only 
cffr^rei which has small annihilation cross section. Thus we can see that the appropriate contribution 
comes from the tt-channel and /i/i-channel. The red, green and blue points correspond to each range 
of sina as written in Fig. 0] In the case of (ys)n = 0.5, the mass relation 2Mi m m# is required 
since the coupling (ys)n is not so large and the annihilation cross section al s v Te \ has to be enhanced 
due to a resonance. The resonance relation is not required for the right side one which corresponds 
to (vs)n — 1-0 if sina is large since the /i/z-channel is effective in this case. There is no allowed 
parameter space in Mi < m h region because the main channel is z/ 1± ^ 1± — > bb and the contribution 
to the annihilation cross section is too small. 



4 Direct Detection 

The DM candidates i>i± interact with quarks via Higgs exchange. Thus it is possible to explore the 
DM in direct detection experiments like XENON100 [33]. The Spin Independent (SI) elastic cross 
section ogi with nucleon N is given by 

a N _ /4m ( 1 1 V Uvs)u musing cos a ^^n\ ^ ^ 
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Figure 5: The comparison with XENON100 experiment. The left figure is for (ys) n = 0.5 and the 
right one is for (ys)u — 1-0. The parameter choice is same as Fig. HJ 

where //dm = (Mf 1 + m^ 1 ) 1 is the DM-nucleon reduced mass and the heavy Higgs contribution 
is neglected. The parameters which imply the contribution of each quark to nucleon mass are 
calculated by the lattice simulation [4*3~||4*4"] as 

f* = 0.023, = 0.032, ff = 0.020, (4.2) 
ft = 0.017, f2 = 0.041, f r ; = 0.020, (4.3) 

for the light quarks and fq = 2/21 (\ — J2 q <3 fq^j f° r ^ ne heavy quarks Q where q < 3 implies the 
summation of the light quarks. The recent another calculation is performed in Ref. |45] . 

The comparison with XENON100 upper bound is shown in Fig. where the other parameters 
are fixed as same as Fig. H] and these correspond to red, green and blue points. The violet dotted 
line is XENON100 (2011) upper bound and the light blue dashed line is XENON100 expected one 
in 2012. We can see that from the figure, the XENON100 (2011) limit excludes M x < 800 GeV in 
the large sin a region for (ys)n — 1-0. The almost excluded region of rather small Mi implies that 
the parameter region of the /i/i-channel is the most effective for the DM annihilation (Fig. H]). The 
other certain region will be verified by the future XENON experiment. In the case of taking into 
account the lightest right-handed neutrino as DM in the original radiative neutrino mass model [B], 
the elastic cross section with nuclei is not obtained at tree level because of the leptophilic feature of 
the DM. However interactions with quarks via Higgses are obtained in the radiative inverse seesaw 
model and hence verification by direct detection of DM is possible. 
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5 Conclusions 



We constructed a radiative inverse seesaw model which generates neutrino masses and includes DMs 
simultaneously, and studied the mixing of the lepton sector and the DM features. The neutrino 
mass matrix was radiatively generated via the inverse seesaw framework. Considering the latest data 
of non-zero #13, we applied the charged lepton mixing effect with almost Maximal Dirac CP phase 
suggested by S. King. We can obtain the neutrino Dirac Yukawa matrix which produces the best fit 
value of 812 on the diagonal basis of the right-handed neutrinos and additional fermions and non-zero 
#13 comes from the charged lepton mixing matrix. The size of the neutrino Yukawa couplings is 
severely constrained by LFV at the same time. As a result, the annihilation cross section which 
comes from the Yukawa interaction becomes ineffective, however we found that new interactions via 
Higgs bosons which are independent on LFV. Thus the DM can have the certain annihilation cross 
section due to the interactions with Higgses. Verification of the model is also possible by direct 
detection of DM through the interaction with Higgses. In particular, the region of the large mixing 
sin a will be testable by the exposure of the XENON100 (2012) experiment. Therefore it is favored 
to be Higgs portal DM in the radiative inverse seesaw model from the view point of avoiding the 
LFV constraint and obtaining the proper detection rate by direct search of DM. 
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